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Abstract 

We consider the large deviation function for a classical harmonic chain composed of N particles 
driven at the end points by heat reservoirs, first derived in the quantum regime by Saito and Dhar 



and Kundu et al. [3|. Within a Langevin 



l|] and in the classical regime by Saito and Dhar 
description we perform this calculation on the basis of a standard path integral calculation in 
Fourier space. The cumulant generating function yielding the large deviation function is given in 
terms of a transmission Green's function and is consistent with the fluctuation theorem. We find a 
simple expression for the tails of the heat distribution which turn out to decay exponentially. We, 
moreover, consider an extension of a single particle model suggested by Derrida and Brunet [J| and 
discuss the two-particle case. We also discuss the limit for large and present a closed expression 
for the cumulant generating function. Finally, we present a derivation of the fluctuation theorem 
on the basis of a Fokker-Planck description. This result is not restricted to the harmonic case but 
is valid for a general interaction potential between the particles. 

PACS numbers: 05.40.-a, 05.70.Ln 
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I. INTRODUCTION 



There is a current interest in the thermodynamics and statistical mechanics of fluctuating 
systems in contact with heat reservoirs and driven by external forces. The current focus 
stems from the recent possibility of direct manipulation of nano-systems and bio-molecules. 
These techniques permit direct experimental access to the probab ility distribution functions 
for the work or for the heat exchanged with the environment js-lSj. These methods have 
also yielded access to the experimental verification of the recent fluctuation theorems which 
relate the probability of observing entropy- generated trajectories with that of observing 



entropy-consuming trajectories 



16 



33|. 



In recent works we studied the motion of a Brownian particle in a general potential with 
a view to the distribution function for the heat exchange with the surroundings 3J] and a 



single bound Brownian particle driven by two heat reservoirs [35|. In the present paper we 



consider the harmonic chain driven by heat reservoirs at temperatures Ti and Tj 



N 



36 



42- 



49j. Here the distribution of positions and momenta is given by a Gaussian form with a 



correlation matrix with elements given by the static position and momentum correlations 



36|. 



Owing to the current interest in fluctuation theorems the linear chain has recently been 
addressed again by Saito and Dhar ^ and by Kundu et al. ; see also jl| for a treatment 
in the quantum regime. Using a path integral formulation, Fourier series, and analyzing the 
resulting energy transmission matrices, these authors derive an expression for the cumulant 
generating function, in the following denoted CGF, for the heat transfer in terms of a 
transmission Green's function T{u). The large deviation function for the heat transfer, 
denoted LDF, then follows by a Legendre transformation of the CGF; _for definitions, see 
later. The expression is in accordance with the fiuctuation theorem fl7 . 



19 



21 



23I-26I. 



In the present paper we consider four issues: i) the CGF for the harmonic chain, ii) the 
LDF for the chain and the exponential tails in the heat distribution, iii) the CGF for an 
extension of a model by Derrida and Brunet ^, iv) the CGF for the harmonic chain in the 
large limit, where is the number of particles, and v) a derivation of the fiuctuation 
theorem on the basis of a Fokker-Planck description. 

For the purpose of the analysis in ii) - iv) we have within a Langevin scheme performed a 
calculation of the CGF, including explicit expressions for the transmission Green function. 
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At the technical level we, moreover, unlike Kundu et al. jsl, make use of Fourier transforms 
throughout the calculation and diagonalize explicitly T{uj) expressing the CGF in terms of 
the eigenvalues. For the benefit of the reader and the continuity of the paper we have chosen 
to include this analysis in the main part of the paper. 

We discuss the tails in the heat distribution and exemplify this feature both for the 
extended Derrida-Brunet model and the particle chain. We consider the CGF and LDF 



for an extension of a single particle model suggested by Derrida and Brunet j4| and the CGF 
in the two-particle case. We, moreover, analyze the asymptotic large limit and present a 
closed expression for the CGF. 

Finally, as a related and more formal issue we present a derivation of the fluctuation 
theorem on the basis of a Fokker-Planck description of a chain. As a bonus we are able to 
prove that the fluctuation theorem holds for chains with general interaction potentials and 
with several heath baths at different temperatures. 

For reference we present below the results of Saito and Dhar ^ and Kundu et al. jsj, 
also presented in the present paper. Denoting the model-dependent transmission Greens 
function by T{uj), the CGF /i(A) for the characteristic function for the transferred heat Q{t) 
in the time interval t, is given by the following expressions: 

(exp(Ag(t))) = exp(t/i(A)), (1.1) 

MA) = -^ / + (1-2) 

/(A) = TiTjvA(l/Ti - 1/Tn - A). (1.3) 

Here Ti and T^r denote the reservoir temperatures and the form of /(A) ensures the validity 
of the fluctuation theorem 

/i(A)=/i(l/Ti-l/Tjv-A). (1.4) 

As a new result we present below the CGF in the asymptotic large A^ limit. Here F 
denotes the reservoir damping and k the spring constant. The CGF is given by 

8r/€-i/2 sin(p/2)sin(p)/(A)" 



K^) = - — v/kcos(p/2) In 
Jo 27r 



:i.5) 



1 + A{T^/k) sin2(p/2) 

Finally, we note that the mathematical background for the present paper is provided by 



large deviation theory, see Refs. 



37^411 
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The paper is organized in the following manner. In Sec. |TT] we present the harmonic 
chain. In Sec. IIIII we set up the necessary analysis. In Sec. |IV] we present a derivation of 
the CGF. The general properties of such a function are discussed in Sec. El where we also 
consider the tails of the heat distribution, the specific cases of a bound Brownian particle, 
a two-particle chain, and the large limit. In Sec. I VI I we discuss a generalization of the 
fluctuation theorem. In Sec. IVIII we present a summary and a conclusion. 



II. HARMONIC CHAIN 



The dynamics of a unit mass harmonic chain composed of particles and, moreover, 
attached to a wall or substrate, is governed by the Hamiltonian 

^ = ^ + f - + f + -^)' (2-1) 

n=l n=l 

where u„ and pn denotes displacements and momenta, respectively; k, is the spring constant. 
The equation of motion for the bulk particles and the end particles driven by the heat 
reservoirs at temperatures Ti and T/v with associated damping F are given by 

dur 



r = Pn, (2.2) 
K{Un+l + Un^i - 2n„,), Tl = 2, ■ ■ ■ iV - 1, (2.3) 



dt 
dpf 

~dt 

^ = -Tp, + k{u2 - 2ui) + ^i, (2.4) 

= -TpN + n{uN-i - 2un) + ^Af, (2.5) 

with noise correlations and strengths 

{mUt')) = ^iS{t-t'), (2.6) 
{^NmN{t')) = Ar,6{t-t'), (2.7) 

Ai = 2rri, (2.8) 

Ajv = 2TTn; (2.9) 

in equilibrium Ai = Aat = A and detailed balance implies A = 2rT, where T is the common 
temperature of the reservoirs. 



5 



Focussing on the reservoir at temperature Ti the fluctuating force is given by —Tpi + 
and, correspondingly, the rate of work or heat flux has the form, denoting Q = Qi, 

dQ 



dt 



(2.10) 



The central quantity in the analysis is, however, the total heat transmitted during a flnite 
time interval t, i.e.. 



Q{t) 



drpi(r)(-rpi(r)+6(r)); 



(2.11) 



note that strictly speaking only the time scaled heat Q{t)/t has large deviation properties. 



see e.g. Refs. 



37 



4l|. 



The heat Q{t) is fluctuating and the issue is to determine its probability distribution 
P{Q,t) = {6{Q — Q{t))); here (■ ■ ■) denotes an average with respect to and ^n- In terms 
of the characteristic function {exp{XQ{t)) we have by a Laplace transform 55| 

■ = dX 



PiQ,t) 



27ri 



(2.12) 



The chain attached to the substrate at the ends and driven by heat reservoirs is depicted in 

Fig.m 



III. ANALYSIS 



The heat reservoirs drive the chain into a stationary state. Since the heat is transported 
ballistically the only damping mechanism is associated with the heat reservoirs and the only 
time scale is given by l/F. Consequently, at long times compared with 1/F we can neglect 
the initial preparation of the chain and analyze the problems in terms of Fourier transforms. 
Thus introducing the Fourier transform 

Mt) = [ %-'-'u4u), (3.1) 
J 27r 

the equations of motion (12. 2p to (12. 5p and noise correlations (12. 6p to (12. 7p take the form 

N 

'^nLi^)Umiuj) = + ^n7V^7v(w), (3.2) 

m=l 

(ei(a;)6(w')) = 27rAi5(a; + u'), (3.3) 

{^N{to)iN{to')) = 2nA^5{cu + cu'). (3.4) 
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Here the inverse Green's function G^^ {u) is a symmetrical tridiagonal matrix with elements 



where 



G^lico) = = n, (3.5) 

G-^{u) = n, n = 2,---iV-l, (3.6) 
G~l+^{uj) = G-^_^{uj) = (3.7) 

n = -uj"^ + 2/4 - iTu, (3.8) 
n = -u^ + 2K- (3.9) 

note that for a free chain we have Vl = — + n — iVu. 
Propagating bulk solutions have the form 

Un(ou) = Aexp{ipn) + Bexp{—ipn), (3.10) 
oj^ = 4:Ksm^{p/2), (3.11) 

where p is confined to the first Brillouin zone \p\ < vr and, correspondingly, < 2^/K,. 
Imposing the noisy drive we readily determine the coefficients A and B and infer the solutions 

Un{uj) = Gnl{uj)^l{uj) + G^TV (w)^Ar (w) , Pn{uj) = {-iuj)Un{(^) , (3.12) 

where the Green function components are given by 

nsm{N -n)p- Ksm{N -n-l)p 
<-^nU^^j 'Diuj) ' [o.i-6) 

nsm{n-l)p- Ksm{n-2)p 
D{u) 

D{uj) = n'^ sm{N - l)p - 2Kn sm{N - 2)p + sm{N - 3)p. (3.15) 

The displacement m„ is thus driven by stochastically excited lattice waves (phonons) prop- 
agating towards the site from the end points; D{uj) = yield the damped mode spectrum. 
Also, from the definition of G~^{uj) we deduce the relationship G^{uj) — G^.^{—uj) = 



-2iuT6nm{Sni + ^un) and by multiplication the Schwinger identity [50 1 



Gnmi(^) - Gnmi(^y = 2ia;r [G„i (u;)Gi„ (w) * + GnNi(^)GNm{(^y]- (3.16) 

In the absence of the heat reservoirs energy is conserved, i.e., dH/dt = 0, where H is 
^iven by (12. ip . Coupling the reservoirs to the chain we have dH/dt = dQi/dt + dQ^/dt, 



where dQ^/dt is the heat flux from the reservoir at temperature T^. Averaging we have 
for the mean heat fluxes (dQi/dt) = —{dQN/dt), expressing the energy balance; the mean 
input flux at n = 1 is equal to the mean output flux at n = A^. 

Using fl2.10p . inserting f l3.12p . averaging over the noises fl3.3p and fl3.4p . using the prop- 
erties of the Green's function (13.130 and (13.140 . and the identity (I3.16p . we obtain for the 
mean transferred heat in time t 

{Q{t)) = t(Ai - A^)r / pu;'\G^^{u)\'. (3.17) 

J 27r 

Here the central model dependent quantity is the end-to-end Greens function Gin{u}); we 
note that the mean heat vanishes for Ai = A^. We also note that the transferred mean 
heat rate q = {Q)/t is given by 

g = r(Ti-(p?)); (3.18) 



see Ref. |36|. Here {pf) is the average kinetic temperature of the flrst particle in the steady 
state; note that in equilibrium q = and {pf) = Ti in accordance with the equipartition the- 



orem [51(1 . The relation (I3.18P follows from (I3.17P by inserting the Greens function solution 
of the equations of motion and using the identity (I3.16p . 

The expression (I3.18P also follows from the equivalent Fokker-Planck approach to the 
harmonic chain which we discuss below. Considering the deflnition of the n-th moment of 
the heat transfer in time t 

{Q'\t)) = j dQdudpQ''Piu,p,Q,t), (3.19) 

and referring to (16. 6p in Sec. IVIl the Fokker-Planck equation for the joint distribution 
P{u,p,Q,t) in the case of two reservoirs implies 

^ = Tn{{T, - pDQ--') + Tn{n - l)Ti(p?Q"-2). (3.20) 

These equations of motion are part of a hierarchy relating the n-th moment to correlations 
of the lower moments with (pf) and have to be completed by equations of motions for the 
correlations Without further assumptions this hierarchy will in general not ter- 

minate and simply represents a reformulation. We note, however, that for the flrst moment 
for n = 1 the second term in (I3.20p vanishes and we obtain a closed equation yielding (13.181) . 
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For the fluctuating heat transferred in time t we obtain, using (12. lip and inserting fl3.12p . 
the expression 

Q{t) = I ^^^(^-^') (eiM e^M)M(a;,a;') (^^^^'^^J^^ (3-21) 

The heat transfer is a fluctuating quantity depending bilinearly on the reservoir noises .^i 
and C,N- The matrix elements in the symmetrical form fl3.2ip are given by 

Mn(a;, u') = -rA{u)A{u'y + {l/2){A{u) + A{u'y), (3.22) 

M22{uj,uj') = -TB{uj)B{uj')\ (3.23) 

Mi2(w, io') = -TA{io)B{uj'y + {l/2)B{uj'y, (3.24) 

M2i{uj, Lo') = -TB{uj)A{uj'y + {l/2)B{uj), (3.25) 

where we have introduced the notation 

A{uj) = -ituGuiuj), (3.26) 

B{uj) = -icoGiNiuj); (3.27) 

we note that (I3.16P implies 

A{uj) + A{ujy = 2T[\A{uj)\^ + \B{uj)\% (3.28) 

The dependence on the transfer time t is embodied in the function 

F(c.)=2e--/^^^^^^. (3.29) 

For later purposes we also note that 

F(0) = t, (3.30) 
|F(cj)p = 27rM(cj) for large t (3.31) 

At this stage our calculation differs from Kundu et al. in that we use a Fourier transform 
instead of a Fourier series in the expression (13.211) for the fluctuating heat. The dependence 
on the transfer time t is then incorporated in the function -F(w). 



9 



IV. LARGE DEVIATION FUNCTION 



For large t the mean heat {Q{t)) given by f l3.17p grows hnearly with time. Analyzing 
the higher cumulants {Q{t)"')c, i-e., (Q(tY)r = {QitY) — {Qit))'^, etc., by averaging over the 



noise and applying Wick's theorem 52|, it also follows that they likewise increase linearly 



with time, i.e., (Q(t)'^) r ^ t for large t. We thus infer from the cumulant expansion of the 



characteristic function 



5l|, 



(exp(Ag(t))) = exp 5^ ^(Q(t)")c , (4.1) 
that for large t 

(exp(AQ(t)))=exp(t/i(A)), (4.2) 

where yu(A) is the cumulant generating function, denoted CGF. 

The CGF characterizes the long time heat distribution. From the cumulant expansion 
(14. ip we obtain the relationship 

/d>(A)\ _ mr)c 

Here the definition (14.21) for A = implies 



(4.3) 



/i(O) = 0. (4.4) 

In case the fluctuation theorem is valid we, moreover, have the symmetry 

/i(A)=/i(l/Ti-l/Tjv-A). (4.5) 

Turning to the evaluation of /i(A) we average exp(A(5(t)) with respect to the noises and 
^Tv- In matrix form the Gaussian noise distribution has the form 

m « exp [-\ J ^^e»A-^(a; - c.O^(-a;')) , (4-6) 

where ^{oo) = {^i{uj),^n{^)) and the inverse noise matrix is given by 

A-\u-u')=i^^ ° ]6{u-u'). (4.7) 
V A-i/ 
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Noting from fl3.2ip that Q{t) is bilinear in ^ and using the identities 52 1 

(exp(-(l/2)f50) = det(J + AB)"^/^, (4.8) 

det{A) = exp(Trln(A)), (4.9) 

we obtain for the CGF 

^(A) = -^Trln(/-2AAFM). (4.10) 

In the remaining part of this section the present calculation differs from Kundu et al. 
in that we owing to the nondiagonal character of F must expand fi in order to implement 
the large t limit. Thus expanding the log according to ln(l + x) = Yln=ii~^)"^~^^^^ /''^ 
tracing term by term we have 

\n+l 



Z^^^) = E ^^(-2A)"Tr(AFM)", (4.11) 

n=l 

and the issue is to determine Tr(AFM)" and complete the sum. From fl3.2ip we obtain 

Tr(AFM)" = /" n ^^(^'^ - ^'^+1)'^^ ( n ^^(^A., c^fc+i) ) , (4.12) 

k=\ \k=\ / 

where Wn+i = ^i- Inserting fl3.29p we notice that since X]fc=i('^fc — exponential 

factors in the product of F functions combine yielding a unit factor. We thus only have to 
retain the sine part, i.e., F{(jj) — )■ 2 sin(a;t/2)/a;. Using f l3.30p and (I3.3ip we have for n = 1, 2 

Tr(AFM) =t /— Tr(AM(w,w)), (4.13) 
J 27r 

Tr(AFM)2 = t / — Tr(AM(u;,u;)AM(w,u;)). (4.14) 
J 27r 

For large t the function F{(jj) oscillates rapidly as a function of u and we have approximately 
Wi ~ Ci;2 ■ ■ ■ ~ oJni i.e., the effective integration range in u space is confined to the domain 
uji = UJ2 = ■ ■ ■ = oJn and only one u integration remains. Using J{du/27i)F{u) = 1 
inspection readily yields 

Tr(AFM)" = t /— Tr((AM(w,w)"), (4.15) 
J 27r 



and the CGF takes the form 

^ -i)"+i , r du 

n=l 
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/^(^) = - 2 E J ^Tr((AM(c., c)"). (4.16) 



In order to complete the calculation we diagonalize the two-by-two matrix AM. Denoting 
the eigenvalues by and a+{uj) we have Tr((AM(ci;, cj)'^) = a+(a;)'^ + a-{u)'^ and 

reconstructing the log we obtain for /i(A) 

/i(A) = -lJ ^[ln(l - 2Xa4uj)) + ln(l - 2A«„(a;))], (4.17) 

The eigenvalues a+(u;) and a+(u;) are determined by the condition det(AM — al) = 0, i.e., 

AiMn{u,u) - a{u) AiMi2{u,uj) 

yielding the quadratic equation 



0, 



- a(AiMn + A^Maa) + AiA;v(MnM22 - M12M21) = 0, 

with roots and a_. In particular 

+ = AiMii + A7VM22, 

= AiA^(MnM22 - M12M21). 

Using the identity (13.281) we obtain the reduced expressions 

Mn{co,co) = -M22iuj,uj) = T\B{uj)\^, 

Mu{uj,co) = M2i{co,ujy = -TA{u)B{uy + (l/2)fi(a;)*. 



I.e., 



«+ + «_ = (Al- AAr)r|fi|2, 

= -AiAAr|5|V4, 



and for the CGF 

fi{X) = -l I ^In [l-2A(Ai-A^)r|i?(a;)|2-A2AiA^|i?(u;)|2] 
Finally, inserting (12. 8 p and (12. 9p the CGF can be expressed in the form 

MA) = -^ / ^Ml + 4r^|5(a;)|V(A)], 



where 



/(A) = TiT^A(-A + 1/Ti - 1/T^) 



(4.18) 



(4.19) 

(4.20) 
(4.21) 

(4.22) 
(4.23) 

(4.24) 
(4.25) 

(4.26) 

(4.27) 

(4.28) 
(4.29) 
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This expression is in agreement with Kundu et ah [3|. Here the form of /(A) ensures that 
the fluctuation theorem (14. 5p holds. The deterministic djTiamics of the chain is entirely 
embodied in the momentum Green's function B{u). 



V. DISCUSSION 

Here we discuss four issues: i) the branch cut structure in /i(A) and ensuing exponen- 
tial tails in the heat distribution P{Q/t), ii) a single bound Brownian particle coupled to 
two reservoirs, iii) a two particle chain coupled to heat reservoirs, and iv) an asymptotic 
expression for /i(A) in the large N limit. 

A. Exponential tails 

By inspection of the general expression (14.271) for the CGF we infer that /i(A) has the form 
of a downward convex function passing through the origin /i(0) = due to normalization 
and through /i(l/Ti — 1/Tjv) = owing to the fluctuation theorem. Since the argument in 
the log in fl4.27l) must be positive we infer the condition 

m > -if4i^' ^^-^^ 

I I max 

where |-B|max is the maximum value of in the u range. By means of algebraic and 

trigonometric manipulations it can be shown that |-B(ci;)p is bounded by l/4r^, for details 
see appendix [XJ and consequently, /(A) > —1. By analyzing the expression for /(A) in 
f l4.29p one easily finds that this bound is satisfied for A_ < A < A+, where the branch points 
A± in /i(A) are given by 

A+ = 1/Ti, (5.2) 
A_ = 1/Tn. (5.3) 

In Fig. [2] we have depicted the CGF given by ( Km for the case Ti = 10, Tn = 12, F = 2, 
K = 1, and iV = 10. 

At large times the heat distribution function follows from (I2.12p . i.e., 

fie-AQ,tM(A)^ (5.4) 

-ioo ^TTZ 
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and the rate function or large deviation function F{q) is given by 



P(g) 

Q 



~ e 



-tF{q) 



(5.5) 
(5.6) 



Since fj,(X) is differentiable, strictly convex, and steep at the boundaries the Gartner-Ellis 



theorem 



37 



41| implies that the LDF is given by the Legendre transform 



F{q) = sup^{q\- 



(5.7) 



or 



P(g,t) ~e 



t(MA*)-A*g) 



where A* is determined by 



/i'(A*) = q, 



and we find the LDF 



F(g) = -/i(A*) + Ay(A*). 



For F{q) the fluctuation theorem has the form 



Fiq)-Fi-q) = qil/n-l/T^). 



(51 



(5.9) 



(5.10) 



(5.11) 



Note that the LDF also follows from a heuristic saddle point argument, see 26|. In Fig. [3] 
we have depicted —F{q) for the case Ti = 10, T/v = 12, F = 2, k = 1, and = 10. 
Replacing yu(A) by the parabolic approximation 



/ipar(A) = q\{TiTN\ + Ti - T, 



(5.12) 



where q is given by flS.lSp we obtain for F{q) 



^par(g) 



{q-qf{T,-T, 



N 



(5.13) 



in accordance with (15.111) . For the heat distribution we obtain the displaced Gaussian 
P{q) oc exp(— tFpar(5')); this also follows from general large deviation theory 



37 



41|. 
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Deforming the contour in the integral (15. 4p to pass along the real axis we pick up branch 
cut contributions in /i(A). Heuristically, we conclude that for large |g| 

F{q) ~ A+g, for g > , (5.14) 
F{q) ~ |A_||g|, for g < 0, (5.15) 

where A+ and A-_ have been defined above. This also follows directly form the Legendre 
transformation since fi{X) is defined on a compact support. The linear behavior is confirmed 
by the plot of F{q) for our particular choice of the parameter set, see Fig. |3l The heat 
distribution thus exhibits exponential tails for large |g|, i.e., 

P{q) oc exp(-A+gt) for g > 0, (5.16) 
P(g) oc exp(-|A_||g|t) for g < 0, (5.17) 

with A+ and A„ given by (15. 2 p and (15. 3p . It is interesting that the tails are determined 
only by the reservoir temperatures. Finally, we note that the exponential tails in P{q) also 



follows from large deviation theory since fi is bounded by A±, see Refs. 37|, |41 |. 



B. Bound Brownian particle 



In an interesting paper Derrida and Brunet 4] considered a single Brownian particle 
driven by two reservoirs at distinct temperatures and presented an explicit expression for 
the CGF /i(A). This toy model has also been discussed by Visco |53| who considered next 
leading term and the role of initial conditions; see also Farago 54 1. 



In a previous paper 



35| we considered an extension of this model to the case of a single 



particle attached harmonically to a substrate with spring constant k using the simple method 
devised by Derrida and Brunet. We found that the CGF is independent of k, indicating that 
the deterministic character of the spring does not influence the statistical properties of the 
long time heat transfer. Here we consider as an illustration the same problem within the 
present scheme and recover a CGF independent of k. The configuration is shown in Fig. |H 
Associating the damping constants Fi and r2 with the two reservoirs the equation of 
motion take the form 



du 

dp 
'dt 



-(Ti + T2)p- + + 



(5.18) 
(5.19) 
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with noise correlations 

{mUt')) = 2r^T^6{t-t'), (5.20) 
{Ut)Ut')) = 2r2T26{t-t'). (5.21) 

In Fourier space we obtain the solution 

p{u) = Biu){Uoo) + U^)), (5.22) 

where 
and 

note that B{u) satisfies the Schwinger identity (I3.28p . i.e., 

B{uj) + B{ujy = 2(Ti + T2)\B{u)\\ (5.25) 

The heat flux from the reservoir at temperature Ti is 

^=p{-T,p + p^,), (5.26) 

and we obtain from (15.220 . fl3.2ip . and fl5.25p the diagonal matrix elements 

Mu{u,uj)= T2\B{u)\', (5.27) 

M22{uj,uj) = -Ti\B{uj)\\ (5.28) 

Mi2(u;,a;) = -T^\BiLo)\^ + (l/2)5(u;)*, (5.29) 

M2i{io,uj) = -Ti\B{io)\^ + {l/2)B{uj). (5.30) 

Following the prescription in Sec. [IV]the eigenvalue equation imply 

a+ + a_ = 2rir2(Ti-r2)|S|2, (5.31) 
a+a^ = -TiT2TiT2\B\'^, (5.32) 



and we obtain the CGF 



^^^^ = ~l J ^ln[l + 4rir2|5(a;)|V(A)], (5.33) 
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where is given by fICTD and /(A) by fOOj) . 



A straightforward evaluation of (15.331) using the integral 55 1 



— m 

27r 



(5.34) 



yields 



/i(A) = (l/2)(a+ + a„ -64 



(5.35) 



where 



:i/2)((ri 



(5.36) 



6| = (l/2)(4rir2/(A) + (Fi + 12)' -2k± v/(4rir2/(A) + (Ii + V^f - 2kY - Ak''. (5.37) 

Further inspection shows, however, that the combination a+ + a_ — 6+ — 6_ is independent 
of the spring constant k, as already shown in jssi, and we obtain 



/i(A) = (1/2) [fi + 12 - v/(r7+T^F+4fjV(A) 



(5.38) 



Introducing /(A), as defined in (14.291) . we can express (15.381) in the form 



^(A) 

where the branch points are given by 
1 



v/rir2Tir2v'(A+-A)(A-A_), 



(5.39) 



Ah 



1/Ti - 1/T2 ± v/(i/Ti - i/T2)2 + (Fi + r2)Vrir2TiT2 



(5.40) 



We note that |i?(a;)p given by (I5.24p has a two-peak structure with maximum value l/(ri + 
r2)^ at oj = and that the expressions (I5.33P and (15.401) are in accordance with the 

general properties of /i(A) discussed above. Finally, using (15. 9p and (I5.10p we obtain for the 
large deviation function F{q), 



F{q) = -(1/2) Ti + 12 - g(A+ + A_) - (A+ - A_) v/fjWVf? 



(5.41) 



which yields a heat distribution P{q) in accordance with the general discussion in Sec. |V] 
A with exponential tails in the heat distribution; for more details regarding this model, see 



35| 
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C. Two particle chain 



As an illustration of the general scheme presented here we briefly consider the case of a 
chain composed of two particles; the configuration is depicted in Fig. O Setting = 2 we 
have from ([S3D, dSH, and (KTI} 



\B{u)\' 



|fi(a;)2 - K 



212 ■ 



2k — iTu, 



(5.42) 
(5.43) 



and we obtain from the general expression (14.271) the CGF 



/x(A) 



du 
2^ 



In 



(5.44) 



We have been unable to reduce the expression (I5.44p further but note that for k = the LDF 
/i(A) = for all A, corresponding to two independent equilibrium systems at temperatures Ti 
and T2. We also remark that decoupling the chain from the walls, corresponding to setting 



-OJ 



K-iTu, we obtain AT^K^uy{\n^-K^\^) = AT^Ky{{u^+T^)\u^-2K+tTu\^). In the 
limit of a stiff chain, corresponding to k — 00, we have 4r^K^u;^/(|r2^ — k^P) — F^/ (w^ + F^), 
i.e., the case of a single unbound particle coupled to two reservoirs, see Sec. IVl B. 



D. N particle chain 

In the limit of large we present below an asymptotic expression for the CGF. For 
general A^ we have from (l3ll|) . (IXTID . (I37[5|) . and Km 

1^/ x|2^ 4/t3 sin^(p/2)sin^(p) 

' ^ ^' \Q^sm{N -l)p-2KQsm{N -2)p + n^sm{N -3)p\^' ^' ' 

which expanding the denominator can be expressed in the form 

' ^ ^' L{p) + K{p)cos{2Np- (t){p)y ^ ' 



where 



a = (1 — a^) cos{p) — 2ia, (5-47) 

6= (l + a^)sin(p), (5.48) 

L=|a|2 + |fe|^ (5.49) 
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M=|6|2-|ap, (5.50) 

C = ab* + a*b, (5.51) 

K = VM2 + C2, (5.52) 

a = (2r/v^) sm(p/2), (5.53) 

tan0 = C/M. (5.54) 

By inspection we note that |-B(a;)p displays an oscillatory structure with approximate period 
Tf/N, reflecting the resonance structure of the propagating lattice waves in the chain. The 
oscillations are modulated by the slowly varying functions of p, sin(j») and sin(p/2). Further 
inspection of f l5.46p shows that the maxima are given by 

\B(u)\^ - 8^"'sin^(p/2)sin^(p) 

where a little analysis implies that |-B(w)|max locks onto l/AT^, corroborating the demon- 
stration of the upper bound in appendix |X1 The lower bound of the oscillatory structure is, 
correspondingly, given by the envelope 

iR/,x|2 8/t-^sin^(p/2)sin^(p) 

the structure is for = 10 depicted in Fig. |9l The positions of the maxima and minima are 
given by the implicit conditions 2Np — (f){p) = ir (mod 2tt) and 2Np — (f){p) = (mod 27r), 
respectively. 

In p space, using du = ^/Kcos{p/2)dp, and inserting (15. 46 p . the CGF given by f l4.27p 
takes the form 

ler^K-i sm\p/2) sm\p)f{Xy 



r dp 

f^W = - — v^cos(p/2)ln 
Jo 



1 + 



L{p) + K{p) cos{2Np - (pip)) 



(5.57) 



In the large limit the rapid oscillations in 



period. Using the integral 55|, see also ref. 

dp 1 



B 



47| 



2 allows us to integrate separately over each 



/o 27ra + 6cos(p) a/o^^^^^^ 
we thus obtain the following approximate form of 

8k~"^ sin^(p/2) sin^(p) 



(5.58) 
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Further, inserting L and K from fl5.49p and fl5.52p we obtain 

2 sin(p/2) sin(p) 



\B\ 



, approx ry^l + 4(rVK)sin2(p/2)' 

and for the CGF in the hmit N ^ oo 



(5.60) 



/i(A) 



''dp 
2^ 



K cos(p/2) In 



1 + 



8rK-^/2 sin(p/2) sin(p)/(A) 
1 + 4(12/^) sin2(p/2) 



(5.61) 



15 



2 

max 



i/4r2 



|5 

|2 



2 

cnv' 



In Fig. [6] we have for = 10, F = 2, and k = 1 depicted \B\'^ 

and |-B|approx- We note that |-B|approx smoothly interpolates over the oscillations in \B\^. In 
Fig. [7] we depict /i(A) as a function of /i for N = 2 and for = 10. The other parameters 
are r = 2,K = l, Ti = l, and Tjv = 1. We note the excellent fit already for = 10 and the 
good approximation at small A for N = 2. 

The expression (I5.6ip for /i(A) is manifestly independent of in the large A^ limit. This 
implies according to fl4.3p that the cumulants and in particular the mean current also are 
independent of A^. This signals that Fourier's law is not valid for the harmonic chain, see 
e.g. ref. 36|. We also note that the large A^ limit does not correspond to the continuum 
limit; we just increase the number of particles in the chain keeping the lattice distance fixed. 



VI. GENERALIZED FLUCTUATION THEOREM 

In Sees. |IV]and |V]we demonstrated the validity of the fluctuation theorem by an explicit 
evaluation of the CGF for the harmonic chain driven at the end points by heat reservoirs 
at distinct temperatures and considered, moreover, some special cases. Here we put these 
results in a more general framework by considering the Fokker-Planck equation for the 
characteristic function 

C(A,t) = (e^^(*)). (6.1) 
For long times C{X,t) ~ exp(t/i(A)) and we obtain the differential equation 

— = /i(A)a (6.2) 
Expressing the Fokker-Planck equation for C in the form 

— = mc (6.3) 
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we identify the CGF fi{X) as the maximal eigenvalue of the Fokker-Planck operator L. The 
issue is thus to establish the fluctuation theorem symmetry for the maximal eigenvalue. 

We aim at ageneralization of the fluctuation theorem to the case of many heat reser- 
voirs, see also jj]. For that purpose we consider a setup where each particle in the chain 
couples to its own heat reservoir at temperature T„. The configuration is depicted in Fig. [HI 
Generalizing (12.101) the heat flux to the n-th particle is given by 

dQn 



dt 



where the noise is correlated according to 



{Ut)U{t')) = 25^mTTJ{t-t'). 



(6.4) 



(6.5) 



Since the transfer of heat induces a change in the state of the system we must at the outset 
consider the joint distribution P{u,p,Q,t) = P({m„}, {pn}, {Qn},t). The heat distribution 
is then given by P{Q, t)_= J dUndpnP{u,p, Q, t). 



As discussed in ref. 



Ill ], the Fokker-Planck equation for the joint distribution P{u, p, Q, t) 



is derived by considering the heat Qn{t) as an independent dynamical variable, whose time 
evolution is governed by (16.41) . Noting that the noise appearing in this equation is correlated 
to the noise appearing in the equation of motion for the momenta (12. 2112. 51) one can write 



dP 
'dt 



{^.^}H-rE(^.S+|:(p"^) 



"^TnPr, 



d^P 

dQndpn 



(6.6) 



where the Poisson bracket is given by 



N 



n=l 



dP dH dP dH 



dpn dUn dUn Opn 



(6.7) 



see also ref. 
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All reference to the deterministic dynamics of the chain is embodied in the Poisson 
bracket. The remaining terms in (16. 6 P are associated with the transfer of heat. Setting 
d/dQn = —Xn and d'^/dQ'^ = we obtain for the characteristic function C{u,p, {Xn},t) 
defined by the multiple Laplace transform 



n ^ exp(- ^ \nQn)C{u, p, {\n}, t), 



(6.^ 
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the Fokker- Planck equation (16.31) . where the operator L{X) has the form 
L{X)C = {C,H} 



T„-^ + (1 - 2A„T„) — (p„C) + {XniKTn - 1)pI + \nTn)C 



(6.9) 



In the absence of couphng between the particles, i.e., for a vanishing Poisson bracket, 
{C, H} = 0, C is the characteristic function for the the heat transfers to independent 
particles coupled individually to reservoirs at temperature T„. Subject to the transformation 
C = exp(_Et) exp((yf)\l/, where g = (1/2) J^ni'^n — ^/'^Tn)p\, the Schroedinger-like equation 
= describes independent oscillators with spectrum E = — r(ni+n2 + - ■ ■ni\f),ni = 
0, 1, ■ ■ ■. The maximal eigenvalue is given hj E = corresponding to /i = 0, characteristic of 
an equilibrium configuration. Turning on the interaction between the particles the maximal 
eigenvalue will be shifted to a finite value and we obtain a nonvanishing A-dependent CGF. 
The structure of (16. 9 p also allows a simple derivation of a generalized fluctuation theorem; 



see also ref. The first step is to perform a "rotation" exp{H/Tm) with respect to the 
m-th reservoir in combination with a time reversal operator T and define the transformed 
Fokker-Planck operator 

Z(A) = e^/^'"ri:(A)r"^e-^/^™. (6.10) 

In the next step we compare the operator L with the adjoint operator L*. Using (d'^/dpf^)* = 
d'^/dPn and (dpn/dpn)* = —pnd/dpn and shifting the Laplace variables A„ it turns out that 
L and L* become identical and we have the relationship 

L(A) = L*(A), (6.11) 

where 

A„ + A„ = 1/T„ - (6.12) 

Since L{X) is related to L{X) by a unitary transformation we infer that L{X) and L*{X) 
have identical spectra and in particular identical maximal eigenvalues, i.e., the same large 
deviation function, 

/^({An})=M{An}). (6.13) 

The expression (I6.13P together with (I6.12p represents a generalization of the usual fluctuation 
theorem to many reservoirs. In the case of two reservoirs, setting T„ = Ti, = T/v, and 
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A„ = A we obtain the usual fluctuation theorem (14.51) . We note that the above derivation 
holds for any time reversal invariant Hamiltonian, i.e., for any kind of interaction between 
the particles. In our derivation we have also assumed that the maximal eigenvalue is positive. 

Thus, our proof of the fluctuation theorem is more general than the one given in Ref. 
which is restricted to the harmonic chain only. Furthermore, our approach does not require 
a direct evaluation of the CGF, but is based only on the property of the dynamics, as 
expressed by the evolution operator L{X). The previous proof can be readily extended to 
the 3-D case, as long as L{X) has a form as in (16.90 . 



VII. SUMMARY AND CONCLUSION 



In this paper we have discussed a variety of issues regarding the noise driven harmonic 
chain. In Secs JIIII and II VI we performed a calculation of the CGF, recovering the results of 
Kundu et al. [3[, but adding some more details for the purpose of our analysis. In Sec IVl we 
discussed the exponential tails in the heat distribution, the bound single particle model, and 
the two-particle chain case. It is an interesting feature of the tails that the fall-off rate only 
depends on the noise features, i.e., the reservoir temperatures, and not on the dynamical 
properties of the chain such as the spring constant n. In the large limit we have found an 
analytical form for the CGF which excellently interpolates the exact result. This result is 
independent of signalling that Fourier'r law does not hold. Finally, incorporating some of 
our results we have in Sec. IVII within a Fokker-Planck description presented a generalization 
of fluctuation theorem to several reservoirs which holds for any interaction potential. The 
fluctuation theorem simply emerges from a symmetry hidden in the Fokker- Planck operator 
and is therefore not restricted to the linear chain but also holds for a 3D system. 
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Appendix A: Maxima of |S(cj)|^ 



Here we analyze the modulus squared of the function 



, —tcuKsmp , 



where 



D{uj) = sm{N - l)p - 2Kn sin(iV - 2)p + sm{N - 3)p, (A2) 
= -w^ + 2ft; - iVuj, (A3) 
= 4A;sin2(p/2), (A4) 

and demonstrate that it is bounded from above by l/4r^, i.e.,|i?(a;)p < l/(4r^). 
Breaking up |-D(c(;)| in real and imaginary parts, 

^\D{uj)\ = -4:KT^sm'^{p/2)sm{N -l)p + K^sm{N +l)p, (A5) 
3[L>(w)] = -4fi;3/2rsin(p/2)siniVp, (A6) 

inserting \B{ijj)\'^, expressing u in terms of p, using (lA4p . and substituting f lA5P - (lA6p . we 
rephrase the condition \B{p)\'^ < l/(4r^) as 

g{p) = 16K^r^ sin^(p/2) [sin^ A^p — sin^p] 

+ (-4fi:r2 sin(p/2) sin(A^ - l)p + sm{N + 1)^)^ > 0. (A7) 

Expressing the sin(p/2) in terms of cosp, and rearranging terms, g{p) becomes 

g(p) = k'^{-T^ sin Np + 2T^ sm{N - l)p - sm{N - 2)p + k sm{N + l)pf, (A8) 

which is non negative thus demonstrating our assertion. The values of p for which g{p) = 
correspond to the points of maximum for with the exception of p = 0, vr where 

D{uj) = 0. So, has N-1 maxima where = 1/(412), pig,^ 



[1] K. Saito and A. Dhar, Phys. Rev. Lett. 99, 180601 (2007). 

[2] K. Saito and A. Dhar, Phys. Rev. E 83, 041121 (2011). 

[3] A. Kundu, S. Sabhapandit, and A. Dhar, J. Stat. Mech. p. P03007 (2011). 

24 



[4] B. Derrida and E. Brunet, Einstein aujourd'hui (EDP Sciences, Les Ulis, 2005). 

[5] E. Trepagnier, C. Jarzynski, F. Ritort, G. Crooks, C. Bustamante, and J. Liphardt, Proc. 

Natl. Acad. Sci. USA 101, 15038 (2004). 
[6] D. Collin, F. Ritort, C. Jarzynski, S. B. Smith, I. T. Jr, and C. Bustamante, Nature 437, 231 

(2005) . 

[7] C. Tietz, S. Schuler, T. Speck, U. Seifert, and J. Wrachtrup, Phys. Rev. Lett. 97, 050602 

(2006) . 

[8] V. Blickle, T. Speck, L. Helden, U.Seifert, and C. Bechinger, Phys. Rev. Lett. 96, 070603 
(2006). 

[9] G. Wang, E. Sevick, E. Mittag, D. J. Searles, and D. J. Evans, Phys. Rev. Lett. 89, 050601 
(2002). 

[10] A. Imparato and L. Peliti, J. Stat. Mech p. P03005 (2006). 

[11] A. Imparato, L. Peliti, G. Pesce, G. Rusciano, and A. Sasso, Phys. Rev. E 76, 050101R (2007). 
[12] F. Douarche, S. Joubaud, N. B. Garnier, A. Petrosyan, and S. Ciliberto, Phys. Rev. Lett. 97, 
140603 (2006). 

[13] N. Garnier and S. Ciliberto, Phys. Rev. E 71, 060101 (R) (2007). 

[14] A. Imparato, P. Jop, A. Petrosyan, and S. Ciliberto, J. Stat. Mech p. P10017 (2008). 

[15] A. Imparato, F. Sbrana, and M. Vassalh, Europhys. Lett 82, 58006 (2008). 

[16] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 

[17] J. Kurchan, J. Phys. A 31, 3719 (1998). 

[18] G. Gallavotti, Phys. Rev. Lett. 77, 4334 (1996). 

[19] G. E. Crooks, Phys. Rev. E 60, 2721 (1999). 

[20] G. E. Crooks, Phys. Rev. E 61, 2361 (2000). 

[21] U. Seifert, Phys. Rev. Lett. 95, 040602 (2005). 

[22] U. Seifert, Europhys. Lett 70, 36 (2005). 

[23] D. J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. Rev. Lett. 71, 2401 (1993). 

[24] D. J. Evans and D. J. Searles, Phys. Rev. E 50, 1645 (1994). 

[25] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995). 

[26] J. L. Lebowitz and H. Spohn, J. Stat. Phys. 95, 333 (1999). 

[27] P. Gaspard, J. Stat. Phys. 117, 599 (2004). 

[28] A. Imparato and L. Peliti, Phys. Rev. E 74, 026106 (2006). 

25 



[29] R. van Zon and E. G. D. Cohen, Phys. Rev. Lett. 91, 110601 (2003). 

[30] R. van Zon, S. Ciliberto, and E. G. D. Cohen, Phys. Rev. Lett. 92, 130601 (2004). 

[31] R. van Zon and E. G. D. Cohen, Phys. Rev. 67, 046102 (2003). 

[32] R. van Zon and E. G. D. Cohen, Phys. Rev. E 69, 056121 (2004). 

[33] T. Speck and U. Seifert, Eur. Phys. J. B 43, 521 (2005). 

[34] H. C. Fogedby and A. Imparato, J. Phys. A 42, 475004 (2009). 

[35] H. C. Fogedby and A. Imparato, J. Stat. Mech. p. P05015 (2011). 

[36] Z. Rieder, J. L. Lebowitz, and E. Lieb, J. Math. Phys. 8, 1073 (1967). 

[37] H. Touchette, Phys. Rep. 478, 1 (2009). 

[38] M. L Freidhn and A. D. Wentzel, Random Perturbations of Dynamical Systems (2nd ed. 

Springer, New York, 1998). 
[39] R. S. Elhs, Physica D 133, 106 (1999). 
[40] S. R. S. Varadhan, Annals of Probabihty 36, 397 (2008). 

[41] F. den Hollander, Large Deviations, vol. 14 (American Mathematical Society, Providence, R.I., 
2000). 

[42] A. Casher and J. L. Lebowitz, J. Math. Phys. 12, 1701 (1971). 

[43] H. Nakazawa, Prog. Theor. Phys. 39, 236 (1968). 

[44] H. Nakazawa, Prog. Theo. Phys. (Suppl.) 45, 231 (1970). 

[45] A. Dhar, Phys. Rev. Lett. 86, 5882 (2001). 

[46] A. Dhar and D. Roy, J. Stat. Phys. 125, 805 (2006). 

[47] D. Roy and A. Dhar, J. Stat. Phys. 131, 535 (2008). 

[48] S. Lepri, R. Livi, and A. Politi, Phys. Rep. 377, 1 (2003). 

[49] A. Dhar, Adv. Phys. 57, 457 (2008). 

[50] M. C. Wang and G. E. Uhlenbeck, Rev. Mod. Phys 17, 323 (1945). 

[51] L. E. Reichl, A Modern Course in Statistical Physics (Wiley, New York, 1998). 

[52] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Oxford University Press, 

Oxford, 1989). 
[53] P. Visco, J. Stat. Mech. p. P06006 (2006). 
[54] J. Farago, J. Stat. Phys. 107, 781 (2002). 

[55] L S. Gradshteyn and I. M. Ryzhik, Table of Integrals. Series, and Products (Academic Press, 
New York, 1965). 

26 




wall wall 



FIG. 1: We depict a harmonic chain in contact with heat reservoirs at temperatures Ti and T/v- 
The chain is attached to walls or substrates at the ends. The total heat transmitted to the n = 1 
and n = N particles are denoted Qi and Qn, respectively. The spring constant is denoted k. 




-0.1 -0.05 0.05 0.1 



A 

FIG. 2: Cumulant generating function /x(A), as given by (|4.27p for Ti = 10, T/v = 12, F = 2, k = 1, 
N = 10. Inset: zoom of the plot for small value of A. 



27 







\ 

\ 


/ 

/ 












\ 

V ^ 










yr / 
^ 1 

1 

1 

1 

1 


\ \ 
\ ^ 
\ * 










1 

1 

1 

1 

1 


V 

I 

I 

I 

1 










1 

1 
1 
1 

1 


I 

I 

1 

1 
t 










; 

( 
1 
f 
1 


1 

t 
1 
1 

• 







-60 -40 -20 20 40 60 

q 



FIG. 3: Full line: plot of the large deviation function —F as a function of q, as given by (|5.10p for 
Ti = 10, T/v = 12, r = 2, K = 1, = 10. Dotted line: parabolic approximation, (j5.13p . Dashed 
and dotted-dashed line: Linear regime for \q\ » q, the slopes are — l/Ti and 1/T/v, respectively. 
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substrate 



FIG. 4: We depict a harmonically bound particle interacting with heat reservoirs at temperatures 
Ti and T2. The heat transferred to the particle is denoted Qi and Q2, respectively. The particle 
is attached to a substrate with a harmonic spring with force constant k. 
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wall wall 



FIG. 5: We depict a chain composed of two particles interacting with heat reservoirs at temper- 
atures Ti and T2- The chain is attached to walls or substrates at the ends. The heat transferred 
to the particles is denoted Qi and Q2-, respectively. The particle is attached to a substrate with a 
harmonic spring with force constant n. 
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0.07 




P 

FIG. 6: We depict the squared modulus \B\^ given by (j5.46p as function of p in the range < p < vr 
for = 10, r = 2, and k = 1 (blue). We also show the maximum value = l/4r^ (black) 

given by (|5.55p , the envelope |-B|env (black) given by (|5.56p . and the large approximation 
l-S|approx (red, dashed) given by (|5.59p. 
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FIG. 7: We depict in two plots the LDF /i(A) as a function of fi for N = 2 and N = 10, respectively. 
The parameters are F = 2, k = 1, and Ti = T/v = 1- The blue curve is based on the exact expression 
given by (j5.57p . the red plusses are given by the N = oo expression in (j5.6ip . 



32 




wall wall 



FIG. 8: We depict a harmonic chain where the n-th particle is in contact with a heat reservoir 
at temperatures T„. The chain is attached to wahs or substrates at the ends. The total heat 
transmitted to the n-th particle ise denoted Qn The spring constant is denoted k. 
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FIG. 9: Plot of the squared modulus of the momentum Green's function 5 as a function of p, as 
given by dS]) for P = 2, k = 1, iV = 10. 
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